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ABSTRACT 

Gravity modes in the Sun have been the object of a long and difficult search 
in recent decades. Thanks to the data accumulated with the last generation of 
instruments (BiSON, GONG and three helioseismic instruments aboard SoHO), 
scientists have been able to find signatures of their presence. However, the in- 
dividual detection of such modes remains evasive. In this article, we study the 
signal at 220.7 /xHz which is a peak that is present in most of the helioseismic 
data of the last 10 years. This signal has already been identified as being one 
component of a g-mode candidate detected in the GOLF Doppler velocity sig- 
nal. The nature of this peak is studied in particular using the VIRGO/SPM 
instrument aboard SoHO. First we analyse all the available instrumental data 
of VIRGO and SoHO (housekeeping) to reject any possible instrumental origin. 
No relation was found, implying that the signal has a solar origin. Using Monte 
Carlo simulations, we find, with more than 99% confidence level, that the signal 
found in VIRGO/SPM is very unlikely to be due to pure noise. 

Subject headings: Sun: helioseismology - Sun: Oscillations 



1. Introduction 

Helioseismology has probed the in- 
terior of the Sun over the last three 
decades. Combining the information pro- 
vided by several hundred pressure- driven 
modes (p modes), it has been possible 
to put constraints on our knowledge of 
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the structure and the dynamics of the so- 
lar interior (Christensen-Dalsgaar dl I2OO2I : 



Thompson et all l2003l ). Unfortunately, 
only a small fraction of these modes 
reaches the solar core. Moreover, due to 
the increase in sound-speed velocity with 
depth, these modes give little information 
on the deeper layers as they spend less 
time there than in the convection zone. 
Let us take as an example the internal ro- 
tation rate of the Sun: the rotation profile 
i s very well known in the convective zone 
flXhompson et al.lll996l : ISchou et al.]ll998l : 
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Howe et all 12000; lAntia fc Basiil 12000), 
while the uncertainties grow in the ra- 
diative zone and towards the core of the 



(Jimenez et al. 


1994 




Els worth et al. 


ICouvidat et al. 2003 


: Chaplin et al. 



20041 : iGarcia et al.ll2004l ). In order to put 
new constraints inside the solar core other 
kinds of modes are needed: the gravity (g) 
modes. For example, by measuring just a 
few of such modes, information on the core 
rotation rate can un doubtably be obtained 
( IMathur et al.l 120081 ). whereas the deepest 
layers that could be probed using; only p- 



modes are around 0.2 Rq ( iGarcia et al. 
2008cj V 



Gravity (g) modes are buoyancy- driven 
modes that have the advantage of prop- 
agating inside the entire radiative region. 
However, these waves become evanescent 
in the convective zone and reach the so- 
lar surface with tiny amplitudes prevent- 
ing us to detect them e a sily ( see for ex- 



ample iBelkacem et al. I (120081 ) and ref- 



erences therein). Indeed, several claims 
for g-mode de tections have been made 
in the past (IPelache fc Scherrer 1983 : 



Palle fc Roca-Cortes 19881 : Thomson et al. 
I995I ): however, modern and better data 
sets cannot confirm them. 

In 1995 was launched the Solar and 
Heliospheric Observatory (SoHO) 



one 



of whose scientific objectives was the 
detection a nd characteriza t ion o f grav- 



ity modes (IDomingo et al.l 1 19951 ). Re- 
cently, using data from the Global Os- 
cillation a nd Low Frequency ( GOLF) in- 



strument (iGabriel et al.l 119951 ). the sig- 



nature of the asymptotic properties of 
£=1 dipole g modes has be en measured 



with a hi^h confidence level (IGarcia et al. 



( IGarcia et al.l 120061 ) using photometric 
data from the Variability of solar IRra- 
diance and Gravity Oscillations ( VIRGO) 
experiment ( IFrohlich et al.l 1 19951 ). Even 
if certain con straints can be imp osed on 
the struc ture (lOarcia et al.ll2008bl ) and dy- 



namics ( Garcia et al. 20071 ) of the solar 
core thanks to the study of these asymp- 
totic properties, it is extremely important 
to detect individual g modes. After 10 
years of observations, the level of noise 
at 200 /iHz has been established at ^ 4.5 
mm/s when individual peaks are looked 
for and at 1.5 mm/ s when this research i s 
done for multiplets ( lElsworth et al.l 120061 ). 
Indeed some peaks and patterns could 
be identified as potential gravity modes 
or mixed modes above the noise level 



(Gabriel et al. 


2002 




Turck-Chieze et al. 


2OO4J; 


Mathur et al. 


200j) but it has been 



impossible to tag them unambiguously 
with the correct ^, m and n. 

In this paper we analyse a peak around 
220.7 /iHz that has been studied several 
times as part of a g-mode candidate us- 
ing different instru ments on boar d SoHO 



(see for example 


Gabriel et al. 


( 


1999): 


Finsterle & Frohlich 


f2001}jTurck-Chieze et al. 


(200i): Mathur et al 


. (2(MI)) but also 



20071 . l2008al ). This signal was also found 



fror a the theoretic a l side (see for exam- 
ple ICox fc Guzik 1 ( 12004 ) and references 
therein). To do so, we start in section 2 
with a brief description of the helioseismic 
instruments used in this work and we anal- 
yse in detail the data of the VIRGO/SPM 
instruments (section 3). In section 4, we 
look for an instrumental origin for this 
peak without success by analysing all the 
housekeeping parameters of the VIRGO 
package as well as the SoHO pointing. 
Once it is stablished that this peak seems 
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to have a solar origin we check for its 
presence in all the other instruments of 
the VIRGO package (section 5) and in 
the velocity instruments GOLF, MDI and 
GONG (section 6). We then finish by dis- 
cussing its possible nature. 

2. Instrumentation and data analy- 
sis 

The data of VIRGO (Variability of IR- 
radiance and Gravity Oscillation), GOLF 
(Gravity Oscillation at Low Frequencies) 
and MDI (Michelson Doppler Imager) on 
board SoHO (Solar and Hehospheric Ob- 
servatory) satellite have been used in this 
research together with the GONG (Global 
Oscillation Network Group) ground-based 
netwok. 

2.1. SoHO/VIRGO 

The VIRGO package was designed to 
study the characteristics of pressure and 
internal gravity modes by observing irra- 
diance and radiance variations, to measure 
the solar total and spectral irradiance and 
to quantify their variability (Frohlich et al. 
1995, 1997). It is composed of three differ- 
ent types of sensors: 

• Two types of absolute radiome- 
ters (one VIRGO/DIARAD and two 
VIRG0/PM06-V) for the measure- 
ments of solar total irradiance and 
its variations with high accuracy 
and precision. The cadences of 
VIRGO/DIARAD and VIRG0/PM06- 
V are 180s and 60s respectively. 

• Two 3-channel sunphotometers (SPM), 
one permanently exposed to sun Ight 
and another for backup, set at at 402 



nm (blue), 500 nm (green), and 862 
nm (red) , looking at the Sun as a star 
with a 60 s cadence. The bandwidth 
of the filters is 5 nm. 

• One Luminosity Oscillation Imager 
(VIRGO/LOI) for the measurements 
of the radiance in 12 pixels over the 
solar disc. The filter is at 500 nm 
with a bandwidth of 5 nm.The ca- 
dence is 60s. 

In 1998 June, SoHO was lost for sev- 
eral months, but, after a search campaign, 
was finally found and resumed operations 
around 1998 October. The VIRGO data 
after SoHO^s "vacations" show the same 
high quality as before the temporary loss 
of the probe. 

2.2. SoHO/GOLF, SoHO/SOI/MDI 
and GONG 

• SoHO/GOLF is a resonance scatter- 
ing spectrophotometer (Gabriel et al. 
1995, 1997) that measures the line- 
of-sight velocity using the sodium 
doublet, similar to the IRIS and Bi- 
SON ground-based networks. The 
GOLF window was opened in 1996 
January and became fully operative 
by the end of that month. Over the 
following months, occasional mal- 
functions in its rotating polarizing 
elements were noticed that led to the 
decision to stop them in a predeter- 
mined position; truly non-stop ob- 
servations began by 1996 mid- April. 
Since then, GOLF has been contin- 
uously and satisfactorily operating 
in a mode unforeseen before launch, 
showing fewer limitations than antic- 
ipated. The signal, then, consists of 
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two close monochromatic photomet- 
ric measurements in a very narrow 
band (25 mA) on a single wing of the 



calibrated into velocity fUlrich et al. 


2000; 


Garcia et al. 


2005) and is in- 



deed similar in nature to other known 
velocity measurements, such as those 
of IRIS and BiSON (Palle et al. 
1999). The sampling of the GOLF 
data used in this paper is 60 s. Be- 
fore SoHO^s vacations (1998 June), 
GOLF data were obtained in the blue 
wing of the sodium line; thus, after 
the SoHO vacations the GOLF team 
decided change to the red wing of the 
sodium line (see Garcia et al. 2005, 
for the latest report on the GOLF 
instrument). 

• The Solar Oscillations Investiga- 
tion (SOI) uses a Michelson Doppler 
Imager (MDI) type of instrument 
( IScherrer et al.lll995l ). MDI consists 
of a pair of tunable Michelson in- 
terferometers, which image the Sun 
onto a 1024 pixel x 1024 pixel CCD 
camera in five wavelengths across the 
Ni i 676.8 nm line. These resolved 
data can be processed by forming 
a weighted combination of the pixel 
signals to yield a p roxy for a Sun-a s- 
a-star response (see Henney ( 19991 )). 



• The ground-based Global Oscilla- 
tion Network Group (GONG; Har- 
vey et al. 1996) consists of six sites, 
with instruments that use the Fourier 
tachometer approach to observe the 
Doppler shift, in the Ni i line, with 
1024 pixel resolution. Here, a Sun- 
as-a-star proxy was formed from a 



simple integration over all pixels with 
a cadence of 60s. 

3. The 220.7 /xHz peak seen in 
VIRGO/SPM 

A long time series of 4098 days of 
VIRGO/SPM data has been used in this 
work starting on 1996 April 11. As our 
purpose is to study the time evolution of 
signals at low frequency, a total of five in- 
dependent subseries of 800 days have been 
computed. Data were available to allow a 
50 day shift up to June 2007. These 66 
overlapped series were used only for plots 
while non-overlapping data were used as 
input for the statistical tests we have car- 
ried on. 

The slow trends in the time series, due 
to the degradation in the instruments and 
long-term solar variability, have been re- 
moved by applying a running mean filter 
of one day. To check whether this filter 
could aflFect the detected signal, we have 
also used a backwards diflFerence filter, in 
which every measured point is substituted 
by the difference of two consecutive points 
= /n+i — fn- To rccovcr the correct 
amplitudes in the power spectrum, this lat- 
ter should be divided by the transfer func- 
tion of the filter Qf^)? defined as follows 
(jGarciafc Ballot I boosi ): 



Q{u) = [2sin{7nyAt)f 



(1) 



u being the frequency and At the sampling 
of the data. 

Both filters gave the same results and, in 
the rest of this paper, we work only with 
the time series filtered by the 1-day run- 
ning mean. 

We therefore computed several power 
density spectra using a Fast Fourier Trans- 
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form (FFT) algorithm and built the time- 
evolution power diagrams used in this 
work. Each of them has been computed 
from the time series which have been ex- 
tended by four equal time intervals of zero 
signal . This oversampling makes it eas- 
ier to detect the bins in which the power 
is concentrated (j Cabriel et al. 2002 1) . We 
also verified that a sine wave fit (SWF), 
computed in steps of 0.0001 /xHz between 
220.5 and 221 muHz^ yields the same re- 
sults. Therefore, we have used the normal 
zero-padded FFT as it is much faster than 
the SWF. 

The time-evolution power diagrams are 
built as follows: The 66 power spectra of 
the overlapped time series are computed 
and plotted vertically using a colour scale 
for the power. The vertical axis is the 
frequency of the power spectra with the 
colour equivalent to the power as indicated 
on the right-hand side of the diagram and 
the horizontal axis the number of the time 
series from to 65, i.e. the time span corre- 
sponding to the time series. Looking these 
time-evolution power diagrams we know at 
which frequency and for how long a sig- 
nal can have enough power to be observed 
above the noise level. 

In figure [1] the time-evolution of the 
three VIRGO/SPM channels are shown for 
the frequency range 220.5-221.0 /xHz. The 
X-axis spans 11 years of the SoHO mission. 
A clear signal is observed in the blue chan- 
nel (top) at around 220.7 /xHz, which is 
stable in time with power that goes from 
6-7 to 16-17 ppm^ / /LiRz. Around time se- 
ries 60 this signal seems to change its fre- 
quency slightly by around 0.3 /xHz. In the 
green and red channels the same continu- 
ous signal is visible as in the blue one but 



with the expected decrease in power with 
wavelength. It is also important to note 
that in aU these VIRGO/SPM channels a 
second high-amplitude signal is visible at 
^ 220.64 /xHz, parallel to the previous one, 
from time series 20 until the last one but 
with a small gap between time series 44 
and 48. 

3.1. Confidence levels and Monte- 
Carlo simulations 

In the previous section we saw that 
there is a persistent signal around the 
target frequency of 220.7 /xHz. Indeed, 
the VIRGO/SPM blue channel power den- 
sity spectrum of the full length time se- 
ries (see Figure [2]) shows the presence of 
a peak at a precise frequency of 220.667 
/iHz above the 90% confidence level com- 
puted in a 10 /iHz win dow following 



Appourchaux et al. I ( 120001 ). 



Using subseries of 800 days and a fre- 
quency window of 10 /xHz, the power 
level above which an observed peak has 
a 90% pr obability not due to n o ise is 8.87o- 



(e.g. see Appourchaux et al. ( 2000l )). In 



the case of zero-padded data, the points 
are no longer independent and are cor- 
related. Therefore, Monte-Carlo simula- 
tions should be used to derive a correc- 
tion for the above-mentioned confidence 
level. In our case, for a padding fac- 
tor of 5 we have added a correction of 



l n(2.8 ) = 1.03 derived by IGabriel et al. 



(120021 ) to the threshold computed using 
non-zero-padded data. In this conditions 
the 90% confidence level at around 9.9cr. 
Using the VIRGO/SPM blue channel, we 
found that the peak we are studying has a 
maximum power in a range between 8.74 
and 10.4a considering only five indepen- 
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ime series number 



221.0 
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Fig. 1.— VIRGO/SPM time-evolution 
power diagram of channels blue, green and 
red (top to bottom respectively) from April 
1996 to June 2007. A stable signal at 220.7 
/iHz is clearly observed. 

dent realizations of 800 days. In Figure [1] 



220 222 
Frequency (/xHz) 



Fig. 2.^ VIRGO/SPM blue channel 
power spectrum density computed with 
4098-day time series starting on 1996 April 
11. The horizontal dotted line corresponds 
to the 90% confidence level that a peak 
above this line would not be due to noise. 

the 90% limit is obtained at around 14 
pprri^ / jiRz (orange colour in Figure [1]). 
This means that, for example, most sub- 
series between the 28th and the 58th have 
the peak above the 90% confidence level, 
as well as other subseries such as those at 
the very beginning of the time-span. It is 
important to notice that in this case the 
G has been averaged over the 66 time se- 
ries and the value of 14 ppm^/fiYLz is an 
averaged magnitude. 

We are interested in knowing the prob- 
ability of having a signal with the same 
properties to those that we have found in 
the VIRGO/SPM blue channel; i.e. a peak 
that is above the 90% level in the full power 
density spectrum of more than 4098 days, 
and that is also present in the five inde- 
pendent subseries of 800 days with similar 
levels to what we have with this instru- 
ment (i.e. not necessary all above a 90% 
confidence level in these individual small 
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subseries but around that level). This lat- 
est condition would be much more restric- 
tive because it means that the peak should 
maintain a certain coherence during the 
full time-span. A Monte Carlo simulation 
of 1 million iterations has been done by 
simulating Gaussian noise time series of 
4000 days that have been cut into five in- 
tervals of 800 days. To speed up the proce- 
dure we have not computed the full spec- 
trum of the 4000 days but only the aver- 
age of the power density spectrum of the 
five independent realizations of 800 days 
(which have an SNR of ^ 9.2 a in the 
VIRGO/SPM blue channel). Thus, the al- 
gorithm looks first for a signal in the av- 
erage spectrum with ^^0.9 times the level 
found in the VIRGO/SPM blue channel 
(8.3(7) and, if it is found, it looks for the 
presence of that signal in the five subseries 
(again with levels of 0.9 times those of 
VIRGO). Any signal with these proper- 
ties found in the 10 /xHz window will be 
fiagged as a positive identification. The re- 
sults show that, in a window of 10 /xHz, the 
220.7 /iHz signal has a likelihood of 99.8% 
(which is reduced to 91.3% if we only con- 
sider the constraint on the averaged spec- 
trum). We have also checked how the like- 
lihood is degraded when a bigger window 
is considered. Thus, for the 20 and 30 /iHz 
windows we obtain 99.6 and 99.4% respec- 
tively. 

We can conclude that it is extremely dif- 
ficult to find a pure noise signal above the 
90% confidence level after ^4000 days and 
with a coherence with time as found in the 
VIRGO/SPM instruments. 



4. Possible instrumental origin of 
the signal inside VIRGO and 
SoHO 

Once this interesting signal has been de- 
tected in VIRGO/SPM the main question 
is to investigate its origin; in other words, 
determine if it is of solar or instrumental 
origin. In this section, we study all the pos- 
sible non-solar origins of this signal, from 
the orbital and pointing corrections of the 
spacecraft to the housekeeping parameters 
(hereafter HK) of the VIRGO package. 

Periodic manoeuvring of the SoHO 
probe at this frequency (220.7 /iHz, i.e. 
a period around 1.25 hours) due to orbital 
adjustments or pointing corrections could 
modulate the signal of the instruments on 
board as a tracking system can produce 
guided frequencies. On the other hand, 
a temperature variation at this frequency 
could also modulate the observed signal. 
These temperature variations could origi- 
nate in the sensor itself or in other instru- 
ment subsystems. 

For all these parameters we follow the 
same analysis we as for the VIRGO data; 
thus, we build the corresponding time- 
evolution power diagrams and we compare 
them with the VIRGO/SPM ones. If the 
220.7 /iHz signal is produced by the tempo- 
ral variations of some of these parameters, 
the time-evolution power diagrams of both 
VIRGO/SPM and the parameter must be 
highly correlated. 

For this purpose we analyse in the fol- 
lowing subsections orbital and pointing 
corrections of the SoHO spacecraft and the 
different HK parameters of the VIRGO 
package that might modulate the signal. 
It is important to note that some of the 
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HK data are in the scientific telemetry of 
VIRGO and have a cadence of 60s, while 
others are in the HK telemetry and have a 
cadence of 180s. 

4.1. Orbital corrections 

The radial distance is reduced to lAU 
by the usual quadratic law Sq = S - r'^ ^ r 
being the spacecraft-to-Sun distance in as- 
tronomical units. This correction normal- 
izes the spectral irradiance to the solar con- 
stant definition and removes signal mod- 
ulations due to movements of the Earth, 
Moon and planets in their orbits. 

The observed radiation S of Sl moving 
blackbody source is 



S — Sq 



(1 



(2) 



where Sq is the radiation in motionless con- 
ditions and V is the speed in units of the 
light speed, c. With SoHO velocity being 
a few 10~^ of the speed-of-light one can 
safely omit terms in v'^ and thus approxi- 
mate the reciprocal formula 



S — S( 



{l-vf 



(l-v^) il-2v) 




- So{l + 2*v) 
(3) 

This Doppler correction removes a tiny 
(10~^), slow (Halo orbit period is 6 months) 
modulation of the measured irradiance. 

In this way, the orbital correction ap- 
plied to the three channels of VIRGO/SPM 
is: 

annel 

(4) 

where "radius" is the spacecraft-to-Sun 
distance in astronomical units and "vel" 
is the radial velocity in units of the speed- 
of-light. 
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Fig. 3. — Time-evolution power diagram 
of the orbital corrections applied to the 
VIRGO data. In addition to the order of 
magnitude which is 6 orders of magnitude 
smaller than VIRGO/SPM no correlation 
is found. 

The orbital parameters (radius and vel) 
are provided by NASA in a 10 minute ca- 
dence and are linearly interpolated to get 
the same 60s as VIRGO/SPM. The time 
series of the orbital correction applied, i.e. 
radius^ {l-\-2'Vel)^ has been analysed in the 
same way as VIRGO/SPM and the result- 
ing time-evolution power diagram diagram 
is shown in Figure [31 The orbital correc- 
tion signal is around 6 orders of magnitude 
smaller than the VIRGO/SPM one and no 
correlation has been found. 

4.2. Spacecraft Pointing 

The three critical fiight dynamics pa- 
rameters are rotations in three dimensions 
around the vehicle's coordinate-system ori- 
gin, the centre of mass. These angles are 
pitch, roll and yaw. Pitch is the rota- 
tion around the lateral or transverse axis. 
Therefore, movements of the spacecraft to 
the north or south of the Sun. Yaw is 
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Fig. 4. — Power time-evolution diagrams 
of the pitch, yaw and roU angles (respec- 
tively from top to bottom) of the SoHO 
spacecraft. Only data up to September 
2002 (time series 32) are available. No 
correlation with the first 32 time series of 
VIRGO/SPM is found. 



the rotation about the vertical axis; thus, 
movements of the spacecraft to the west or 
east of the Sun and, finally, roll is a rota- 
tion around the longitudinal axis i.e. move- 
ments of the spacecraft from the north or 
south to the west or east of the Sun. 

For an instrument that looks at the 
Sun as a star (integrated light) the most 
plain pointing correction would be divided 
by cos{A^yaw'^ + pitch^)^ i.e. the cosine 
of the angle between instrument optical 
axis and the line-of-sight direction. Nev- 
ertheless, this correction was never ap- 
plied to VIRGO/SPM because the correc- 
tion would have been negligible. Also, in 
December 2001, NASA discontinued the 
CDHF (Central Data Handling Facility), 
which was the facility in charge of process- 
ing, producing and distributing the SoHO 
telemetry and the ancillary data products. 
The production of all these data were con- 
tinued in others ways but the production 
of attitude data was stopped in Septem- 
ber 2002. Indeed, when SoHO is in normal 
mode, the attitude follows nominal atti- 
tude well enough for most purposes, and 
because the roll determination had large 
errors (because of certain procedural prob- 
lems) . 

Even knowing that, it would be very un- 
likely that the pointing maneuvres could 
modulate any signal in the SoHO instru- 
ments; the three angles have been analysed 
and their time-evolution diagrams com- 
puted (see Figure [31). 

The available attitude data concerning 
pitch, yaw and roll angles were obtained 
from the NASA archive from 1996 April 
11 to 2002 September 22 and we built 
the time series of the three angles. These 
data have a cadence of 10 minutes and 
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have been used with this samphng rate be- 
cause it is good enough for our purposes. 
The length of these time series enables us 
to get 32 time series of 800 days (each 
shifted 50 days with respect to the previ- 
ous one) . This is approximately half of the 
time-evolution power diagrams used in the 
VIRGO/SPM. This length is sufficient to 
see if any correlation exists between point- 
ing and SPM signals during the common 
period (around 6 years). 

The pitch angle has a very constant 
value of around -3.3 arcmin during the 
time span, with spikes of 5 arcmin and 
only a few of them with higher values, 
between 6 and 13 arcmin. These lat- 
ter are probably due to spacecraft maneu- 
vres. The associated time-evolution dia- 
gram is shown in Figure [4] (top). Some 
power density has been found at a level of 
10~^ {arcminY / /LiHz with no visible corre- 
lation with the VIRGO/SPM signal. 

The yaw angle is zero during practi- 
cally the whole time-span considered with 
some spikes around 1.7 arcmin and only a 
few between 6.8 to 12 arcmin. This yields 
a pure noise time-evolution diagram (see 
Figured] (medium)) with a power density 
of around 10~^ (arcmin)'^ / /xHz with also no 
visible correlation with the VIRGO/SPM 
signal. 

Finally, as we have already said, the 
roll angle does not affect the data achieved 
by instruments that observe the Sun as 
a star (integrated light) but, in any case, 
it has also been analysed. The roll an- 
gle changes following the Earth orbit be- 
tween 7.16 and -7.16 degrees with some 
large spikes that have been removed (the 
roll angle sometimes has large errors) and 
the time-evolution diagram is shown in 



Figured] {bottom). The power density is 
around l(arcmm)^//xHz and, once again, 
no correlation with the measurements of 
VIRGO/SPM has been found. 

4.3. VIRGO/SPM temperatures 

• VIRGO/SPM sensor temperatures 

The most important VIRGO/SPM 
temperature is the temperature sen- 
sor. Each of the three VIRGO/SPM 
channels (blue, green and red) are 
corrected by a quantity proportional 
to each of the temperature sensors 
(sensor blue, green and red). This 
correction is applied in the level 1 
software, so the data we are handling 
are already multiplied by this quan- 
tity. This correction is: 

S PMchannel — channel {T S channel ~'^^^-^^)) 

(5) 

where "channel" means blue, green 
or red; C channel is a constant for each 
channel and TSchannei is the temper- 
ature of each of the three sensors. 

In Figure [5] the time-evolution power 
diagram of the temperature sensor 
of the blue channel is shown (top 
of the right column) . The fluctu- 
ation of this temperature is two or- 
ders of magnitudes smaller that the 
VIRGO/SPM signals at frequencies 
around 220/xHz and no clear corre- 
lation with VIRGO/SPM signals is 
visible. 

• VIRGO/SPM electronic tempera- 
ture 

The temperature of the SPM elec- 
tronics has been also analysed to see 
if there exists some modulation that 
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Fig. 5. — Time-evolution power diagram of the HK temperatures of the VIRGO/SPM pack- 
age. To simphfy the comparison we have repeated, on the same scale, the time-evolution 
diagram of the VIRGO/SPM Blue channel. 



could produce a periodic variation 
in the output voltage of the low- 
noise electrometer amplifiers (or in 
the input current). If this exists, a 
modulation would go to the Voltage 
Frequency Converters (VFC) of the 



Data Acquisition System (DAS) and 
could produce a modulation in the 
output signal. 

The time-evolution power diagram 
for the VIRGO/SPM electronic tem- 
perature is shown in Figure [5] (mid- 
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die of the left column). The power 
density is of the same order as in the 
VIRGO/SPM channels but no corre- 
lation with the signal at 220.7 /xHz 
has been found. 

• Data Acquisition System (DAS) tem- 
perature 

The Data Acquisition System (DAS) 
of VIRGO comprises the Onboard 
Data Handling System (interface for 
telemetry, telecommands and timing 
signals), multiplexers. Voltage Fre- 
quency Converters (VFC) and coun- 
ters. If the DAS temperature, the 
VFC or the counters have a peri- 
odic behaviour, the output number 
of counts could contain that period- 
icity. The DAS temperature time- 
evolution diagram is shown in Fig- 
ure [5] (middle of right column) . The 
power density is an order of magni- 
tude higher than the SPM signal but 
again there is no correlation with the 
signal at 220.7 /xHz . 

• VIRGO/SPM Heatsink and DC/DC 
temperatures 

The temperature variations of the 
VIRGO Heatsink and the VIRGO 
Power Supply (DC/DC) have been 
also analysed for security. The 
Heatsink time-evolution diagram (Fig- 
ure ED (bottom of the left column) is 
ten times smaller than the VIRGO/SPM 
and that corresponding to the DC/DC 
is ten times larger (Figure [5] , bottom 
of the right column). In both cases 
no correlation is found with the 220.7 
/iHz signal. 

Finally, if a signal is the result of 
a certain temperature modulation. 



the temperature variation would be 
higher just where the power of the 
signal is higher. In Figure [6] the 
SPM/Blue and VIRGO HK time- 
evolution power diagrams are plot- 
ted together but only between time 
series 35 to 55, in which the power 
of the 220.7 /xHz signal is stronger in 
the VIRGO/SPM data. This zoom 
helps us to see the darker parts of 
some HK time-evolution diagrams 
produced by the colour scales. None 
of the temperatures analysed in this 
section can explain the observed sig- 
nal at 220.7 /iHz. 

5. The 220.7 //Hz signal in the oth- 
ers VIRGO instruments 

As was mentioned in section 2.1 the 
VIRGO package comprises the SPM Sun- 
photometers and also two types of abso- 
lute radiometers (VIRGO/DIARAD and 
VIRG0/PM06-V) and one Luminosity 
Oscillation Imager (VIRGO/LOI). In this 
section we study the 220.7 /iHz signal in 
these instruments. 

• Luminosity Oscillation Imager (VIRGO/LOI) 

VIRGO/LOI measures the radiance 
in 12 pixels over the solar disc. We 
convert these 12 pixels into one by 
simply adding all of them. From 
the raw time series the same analysis 
as in VIRGO/SPM has been carried 
out. Figure [3 (Top) shows its time- 
evolution diagram. It looks similar to 
the VIRGO/SPM and with the same 
visible signal at 220.7 /xHz. In the 
VIRGO/LOI observations, the signal 
is weaker than in VIRGO/SPM but 
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Fig. 6. — Zoom of the time-evolution power diagrams of the VIRGO/SPM Blue and the 
VIRGO housekeeping analysed in this researtg^ for the time series 35 to 55 where the 220.7 
/xHz has higher amplitudes. This zoom helps to clarify the darker parts of some time- 
evolution diagrams produced by the colour scales. None of these temperatures can explain 
the observed signal at 220.7 /xHz . 
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with the same characteristics, for ex- 
ample, at time series 60 the signal 
slightly changes its frequency. How- 
ever, with this instrument, the peak 
seems to be like a doublet instead of 
only one signal concentrated in a cou- 
ple of bins. 

• Absolute radiometers (VIRGO/DIARAD 
and VIRG0/PM06-V) 

Absolute radiometers use a quite dif- 
ferent technique from that of VIRGO/SPM 
and VIRGO/LOI, which are sili- 
con detectors measuring the spec- 
tral irradiance and the radiance re- 
spectively. Absolute radiometers 
are based on the measurements of 
the heat flux by using an electri- 
cally calibrated heat flux transducer 
to measure the total solar irradi- 
ance (solar constant). Once again, 
from the raw time series we have 
performed the same analysis. Fig- 
ure [71 {medium and bottom) shows 
the results for VIRGO/DIARAD and 
VIRGO /PM06- V. Time-evolution 
power diagrams are similar to the 
previous ones, although the 220.7 
/xHz is weaker in both radiometers 
but the 220.7 /iHz signal is still 
present. 

Analysis using Doppler velocity 
instrumentation 
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Up to now we have not found any instru- 
mental origin for the 220.7 /xHz signal ob- 
served in all the VIRGO package. We can 
now study this region in other helioseismic 
instruments. We will start by analysing 
the signal of the other two instruments on 
board SoHO and we finish by using the 
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Fig. 7. — Time-evolution power dia- 
gram of the Luminosity Oscillation Imager 
(VIRGO/LOI), the VIRGO/DIARAD ab- 
solute radiometer and the VIRG0/PM06- 
V absolute radiometer, respectively from 
top to bottom. All these instruments are 
part of the VIRGO package. 
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GONG ground-based network. 



• GOLF is the other Sun-as-a-star in- 
strument on board SoHO. We have 
analysed the velocity time series 
following the same procedure em- 
ployed in the VIRGO analysis and 
we have computed the time-evolution 
power diagram shown in Figure M 
(top). As mentioned in the intro- 
duction, the 220.7 /xHz signal was 
first observed by GOLF during the 
first years of the mission and it was 
fiagged as a "g-mod e " can didate by 



Turck-Chieze et al. I (120041 ) and, af- 



ter 4182 days, it is still visible as part 
of a quadrupl et above a 90% con - 
fidence level ( IGarcia et al.l l2008al ). 
Figure M ( Top ) shows that the evo- 
lution with time of the signal, al- 
though weaker than in VIRGO, is 
still there. The signal in GOLF has 
an interval between time series 12 
and 19, where it disappears, and it 
corresponds to the place where the 
signal in VIRGO/SPM is the weak- 
est (see figured]). Therefore, we can 
conclude that the 220.7 /xHz is also 
observed in velocity measurements 
using GOLF data but with a smaller 
signal-to-noise ratio (SNR). 

• Disc-averaged MDI velocity signals 
from the calibrated level- 1.4 MDI 
LOI-proxy Doppler images were ob- 
tained using integrated spatially 



weighted masks following iHenney 
fil999l These time series from 
1996 May 25 tiU 2007 October 28 
have been analysed and the time- 
evolution diagram plotted in Figure [8] 
(medium). There are no fingerprints 
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Fig. 8. — Time-evolution power diagram of 
the GOLF, MDI and GONG instruments 
respectively from top to bottom. 
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of the presence of the 220.7 /xHz sig- 
nal in this data set. This could be 
due to a lower SNR in the MDI LOI- 
proxy as compa r ed to GOLF. Indeed, 



Henney et al. I ( 119991 ) showed that 
for the lowest measurable p-modes, 
the GOLF instrument has a higher 
SNR than this particular MDI mask. 
To go further, we need to apply our 
methodology to specifically designed 
g-r node masks such as th ose derived 



bv lWatcher et al. I fl2003l ). 



• The radial velocity of GONG used 
in this work started on 1995 May 7 
and finished on 2006 March 9. These 
series are shifted by a year com- 
pared to the SoHO data but they 
were the longest we could use. The 
disc-integrated data provided by the 
GONG Team are very noisy at low 
frequency and were not suited for our 
studies. Therefore, we have decided 
to use the ^=2 spherical-harmonic se- 
ries. These are optimized for acous- 
tic modes of this degree and they 
have the advantage of having a much 
stabler behaviour at low frequency. 
We preferred these series to the £=1 
mode because the closest theoreti- 
cal frequency to the target frequency 
of 220.7 /xHz corresponds to an £=2 
g mode. In Figure [8] (bottom) the 
time-evolution diagram of GONG is 
shown. Although the SNR at this 
frequencies in the GONG data is very 
low, it seems to be a trace of the 
220.7 /xHz signal in the GONG data, 
especially between series 36 and 52, 
which corresponds to a maximum in 
the GOLF between series 43 and 60 
(7 subseries shift, i.e. ^ 350 days). 



However, looking only at this time- 
evolution diagram, it is impossible 
to disentangle the feature at 220.7 
/xHz from others visible in the anal- 
ysed region. Nevertheless, It would 
be extremely important to be able 
to confirm the visibility of such a 
peak using ground-based data be- 
cause that would directly mean that 
the nature of this 220.7 /iHz signal 
has a solar origin. 

To compare the averaged behavior of 
the 220.7 //Hz signal in the GOLF, GONG 
and VIRGO/SPM data sets we computed 
the coUapsograms of the time-evolution 
power diagram, i.e. to average the 66 power 
spectra used to produce the time-evolution 
power diagrams. The resultant graphs are 
plotted in Figured A similar structure ap- 
pears around the target frequency of 220.7 
/iHz, this peak being the highest in the 
three instruments. However, in the case 
of GONG, it is at noise level. 

7. Conclusions 

In the present paper we have studied a 
peak that appears around the frequency 
of 220.7 /xHz in the VIRGO/SPM data. 
This peak has a more than 90% confidence 
level of not being due to noise in the full 
spectrum of 4098 days. A detailed study 
of its nature revealed that this peak ex- 
isted since the very begining of the mis- 
sion in a continuous way for the last 11 
years and only at the very end of the time 
series considered does it seem to change 
slightly in frequency. By Monte Carlo sim- 
ulations we have computed the confidence 
level of such kinds of behaviour and we 
found that it is really unlikely (more than 
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Fig. 9. — CoUapsograms of the time- 
evolution power diagrams of VIRGO/blue, 
GOLF and GONG. The 220.7 /iHz struc- 
ture is present in the three different in- 
struments although in GONG it is at noise 
level. 

99%) that it is due to a noise with the same 
statistical characteristics as the convective 
noise. Therefore, we checked all the avail- 
able housekeeping data from the VIRGO 
package as well as a detailed analysis of the 
SoHO spacecraft attitude control, looking 
for an instrumental origin. None of these 
studies was able to explain the presence 
of a peak in the region studied. Indeed, 
this study seems to rule out this possi- 
bility. The origin should therefore be so- 
lar. We then studied Doppler velocity data 
from another instrument on board SoHO, 



GOLF, and we found that the peak is also 
present (with lower SNR). Even though 
analysis of data from the GONG ground- 
based network revealed a very noisy spec- 
trum, the highest peak in a 10 /iHz region 
around the 220.7 /iHz signal is precisely 
that peak. However, it is not significant 
enough for us to claim that we have a pos- 
itive detection using this instrument. 

The present study has proved the solar 
origin of the peak at 220.7 /iHz. Two so- 
lar phenomena could be responsible of such 
a peak. The first could be convection, in 
particular granulation motions. However, 
it is very unlikely that a turbulent displace- 
ment of plasma on the solar surface with a 
typical time scale of 10 minutes give a sta- 
ble frequency during more than 10 years in 
the power spectrum of the disc-integrated 
data. On the other hand, gravity modes 
propagate inside the radiative region of the 
Sun and are expected to have long lifetimes 
(at least longer than the period of measure- 
ments). Thus, the properties of the peak 
that we found are similar to those expected 
for a g mode. Using the principle of Ock- 
ham's razor (or the Lex Parsimoniae prin- 
ciple) in which the explanation of any phe- 
nomenon should make as few assumptions 
as possible, we can conclude that if this 
peak is not noise it should be a component 
of a g mode. Analysing in detail the struc- 
ture of this possible g-mode component, 
Figure 2 reveals a peak structure contain- 
ing several bins. Indeed, Figure 1 might 
also show the presence of a parallel com- 
ponent at around 220.64 /iHz with high 
amplitudes in several of the series consid- 
ered. Thus, a possible explanation of such 
behaviour might be the presence of an in- 
ner magnetic field that could slightly split 
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the component of the g-mode multiplet in 
some peaks. Another possibihty might be 
that the g-mode power could be spread into 
several bins as a consequence of a smaller 
than expected lifetime or due to a change 
in the size of the resonant cavity (for exam- 
ple due to a displacement of the position 
of the tachocline during the activity cycle). 
This latter effect is particularly interest- 
ing because it seems that the 220.7 signal 
follows a small change in frequency over 
the entire time-span with the lowest fre- 
quency (220.68 /xHz) reached around time 
series 35 — corresponding to the maximum 
of the activity cycle — and then increas- 
ing the frequency again towards the two 
periods with minimum activity (at the be- 
ginning and the end of the series). In any 
case, assuming a faster rotation in the core 
than in the rest of the radiative envelope 
(as suggested by Garcia et al. 2007), the 
220.7 /xHz peak could be or a component 
of the ^=2, n=-3 g mode, or a component 
of the ^=3, n=-5 or a bitting between this 
latter and the £=5, n=-8. Whatever the 
true answer is, there is still an important 
question to be answered: why is this par- 
ticular peak so excited when there are no 
other visible g-mode components? More 
work will be necessary before solving the 
solar g-mode puzzle. 
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